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Background: The transport through a quantum-scale device may be characterized by the 
^\\ , transmission eigenvalues. These values constitute a junction PIN code where, for example, 

in single-atom metallic contacts the number of transmission channels is also the chemical 

j^ ' valence of the atom. Recently, highly conductive single-molecule junctions (SMJ) with mul- 

^ ! 
I tiplc transport channels have been formed from benzene molecules between Pt electrodes. 

O . Transport through these multi-channel SMJs is a probe of both the bonding properties at 

the lead-molecule interface and of the molecular symmetry. 

g ' Results: Here we utilize a many-body theory that properly describes the complementary 

nature of the charge carrier to calculate transport distributions through Pt-benzene-Pt junc- 
tions. We develop an effective field theory of interacting pi-electrons to accurately model 



■M 



^ . the electrostatic influence of the leads and an ab initio tunneling model to describe the 

lead-molecule bonding. With this state-of-the-art many-body technique we calculate the 
,,_^ ^ transport using the full molecular spectrum and using an 'isolated resonance approximation' 

■"nI" , for the molecular Green's function. 

in 

OO ■ Conclusions We confirm that the number of transmission channels in a SMJ is equal to 

in 

• ' the degeneracy of the relevant molecular orbital. In addition, we demonstrate that the iso- 

^^ I lated resonance approximation is extremely accurate and determine that transport occurs 

predominantly via the HOMO orbital in Pt-benzene-Pt junctions. Finally, we show that 
the transport occurs in a lead-molecule coupling regime where the charge carriers are both 
p\^ • particle-like and wave-like simultaneously, requiring a many-body description. 



INTRODUCTION 

The number of transmission channels for a single-atom contact between two metallic electrodes 
is simply given by the chemical valence of the atom [l|]. Recently 2] it has been determined that 
the number of dominant transmission channels in a single-molecule junction (SMJ) is equal to the 
degeneracy of the molecular orbital [3| closest to the metal Fermi level [4]. In this article, we focus 
on ensembles highly conductive Pt-benzene-Pt junctions [4] in which the lead and molecule are in 
direct contact. 

For a two-terminal single-molecule junction (SMJ), the transmission eigenvalues r„ are eigen- 
values of the elastic transmission matrix [5|] 

T{E)=T^{E)G{E)T^{E)G^{E), (1) 



where G is the retarded Green's function [^ of the SMJ, F" is the tunneling- width matrix describing 
the bonding of the molecule to lead a and the total transmission function T{E) = Tr {T{E)}. The 
number of transmission channels is equal to the rank of the matrix (1), which is in turn limited by 
the ranks of the matrices G and F" 2] ■ The additional two- fold spin degeneracy of each resonance 
is considered implicit throughout this work. 

As indicated by Eq. [H any accurate description of transport requires an accurate description 
of G, which can be calculated using either single-particle or many-body methods. In effective 
single-particle theories, including current implementations of density functional theory (DFT), it 



is often necessary [7l4lOl| to describe the transport problem by considering an "extended molecule," 
composed of the molecule and several electrode atoms. Although this procedure is required in 
order to describe transport at all, it make it difficult, if not impossible, to assign transmission 
eigenchannels to individual molecular resonances since the extended molecule's Green's function 
bears little resemblance to the molecular Green's function. 

We utilize a nonequilibrium many-body theory based on the molecular Dyson equation (MDE) 
[6] to investigate transport distributions of SMJ ensembles. Our MDE theory correctly accounts 
for wave-particle duality of the charge carriers, simultaneously reproducing the key features of both 
the Coulomb blockade and coherent transport regimes, alleviating the necessity of constructing an 
"extended molecule." Consequently, we can unambiguously assign transmission eigenchannels to 
molecular resonances. Conversely, we can also construct a junction's Green's function with only a 
single molecular resonance. The theory and efficacy of this 'isolated resonance approximation' are 
investigated in detail. 



12l ] to transport through SMJs utilized a semi- 



Previous applications of our MDE theory [6|, 111 
empirical Hamiltonian [13| for the vr-electrons, which accurately describes the gas-phase spectra of 
conjugated organic molecules. Although this approach should be adequate to describe molecules 
weakly coupled to metal electrodes, e.g. via thiol linkages, in junctions where the vr-electrons bind 
directly to the metal electrodes [J], the lead-molecule coupling may be so strong that the molecule 
itself is significantly altered, necessitating a more fundamental molecular model. 

To address this issue, we have developed an effective field theory of interacting ir-electrons (vr- 
EFT), in which the form of the molecular Hamiltonian is derived from symmetry principles and 
electromagnetic theory (multipole expansion). The resulting formalism constitutes a state-of-the- 
art many-body theory that provides a realistic description of lead-molecule hybridization and van 
der Waals coupling, as well as the screening of intramolecular interactions by the metal electrodes, 
all of which are essential for a quantitative description of strongly-coupled SMJs \4^. 

The bonding between the tip of electrode a with the molecule is characterized by the tunneling- 
width matrix F", where the rank of F" is equal to the number of covalent bonds formed between the 
two. For example, in a SMJ where a Au electrode bonds to an organic molecule via a thiol group, 
only a single bond is formed, and there is only one non-negligible transmission channel IJ, ll5|. In 
Pt-benzene~Pt junctions, however, each Ft electrode forms multiple bonds to the benzene molecule 
and multiple transmission channels are observed [|4j . In such highly-conductive SMJs the lead and 
molecule are in direct contact and the overlap between the vr-electron system of the molecule and 
all of the atomic- like wavefunctions of the atomically-sharp electrode are relevant. For each Ft tip, 
we consider one s, three p and five d orbitals in our calculations, which represent the evanescent 
tunneling modes in free space outside the apex atom of the tip. This physical model for the leads 
accurately describes the bonding over a wide range of junction configurations. 

In the next section, we outline the relevant aspects of our MDE theory and derive transport 
equations in the isolated resonance approximation. We then outline the details of our atomistic 
lead-molecule coupling approach, in which the electrostatic influence of the leads is treated using 
TT-EFT and the multi-orbital lead-molecule bonding is described using an atomistic model of the 
electrode tip. Finally, the transport distributions for an ensemble of Pt-benzene-Pt junctions are 
shown using both the full molecular Green's function and using the isolated resonance approxima- 
tion. 



MANY-BODY THEORY OF TRANSPORT 

When macroscopic leads are bonded to a single molecule, a SMJ is formed, transforming the 
few-body molecular problem into a full many-body problem. The bare molecular states are dressed 
by interactions with the lead electrons when the SMJ is formed, shifting and broadening them in 

accordance with the lead-molecule coupling. 

n 

Until recently [6] no theory of transport in SMJs was available which properly accounted for the 
particle and wave character of the electron, so that the Coulomb blockade and coherent traiisport 
regimes were considered 'complementary ' [lOJ]. Here, we utilize a many-body MDE theory [g, ll2| 



based on nonequilibrium Green's functions (NEGFs) to investigate transport in multi-channel SMJs 
which correctly accounts for both aspects of the charge carriers. 

In order to calculate transport quantities of interest we must determine the retarded Green's 
function G{E) of the junction, which may be written as 

G{E)=\sE-H^^^,-i:{E)y\ (2) 



(1) ^ f7(2) 



el 



is the molecular Hamiltonian which we formally separate into one-body 



12]. S is an overlap matrix, which in an orthonormal basis reduces to the 



where H^oi=H^^ + H^ 
and two-body terms [6 
identity matrix, and 

^{E) = ^^{E) + ^^{E) + ^c{E), (3) 

is the self-energy, including the effect of both a finite lead-molecule coupling via Srp' and many- 
body interactions via the Coulomb self-energy Y,c{E). The tunneling self-energy matrices are 
related to the tunneling-width matrices by 

r»{E) ^ i (s|(i?) - [s^(i5;)]t) . (4) 

Throughout this work we shall invoke the wide-band limit in which we assume that the tunneling 
widths are energy independent T'^{E) ~ F". 

It is useful to define a molecular Green's function Gnioi(-E') = limpa_j.o+ G{E). In the sequential 
tunneling regime [a], where lead-molecule coherences can be neglected, the molecular Green's 
function within MDE theory is given by 

G^oi(i^) = [^E - Hi^l - J:^S\e)] "' (5) 

where all one-body terms are included in H^^^^ and the Coulomb self-energy S'-''-' accounts for the 
effect of all two-body intramolecular many-body correlations exactly. The full Green's function of 



the SMJ may then be found using the molecular Dyson equation 

G{E) = GrnoiiE) + G^oi{E)AJ:{E)G{E), (6) 

where AS = St + ASc and ASc = Sc — S^ . At room temperature and for smah bias voltages, 



ASc ~ in the cotunnehng regime [6|] (i.e., for nonresonant transport). Furthermore, the inelastic 
transmission probability is negligible compared to the elastic transmission in that limit. 

The molecular Green's function Gmoi is found by exactly diagonalizing the molecular Hamil- 
tonian, including all charge states and excited states of the molecule. Projecting onto a basis of 



relevant atomic orbitals one finds 



1^ 



G^oi(^) - 2^ E_E^,^E^^iQ+ ' (7) 

where V{i') is the probability that the molecular state v is occupied, C{i', v') are many-body matrix 
elements and //moi \y) = Ey\v). In linear-response, V{i')=e~^'^^"~^^"^ / Z, where Z=Y^^ ^-PiE^-N^ti) 
is the grand canonical partition function. 
The rank-1 matrix C(z^, i'') has elements 

where dna annihilates an electron of spin a on the nth atomic orbital of the molecule and i' and v' 
label molecular eigenstates with different charge. The rank of C(z^, z^') in conjunction with Eqs. [6] 
and [7] implies that each molecular resonance v ^ v' contributes at most one transmission channel 
in Eq. [H suggesting that an M-fold degenerate molecular resonance could sustain a maximum of 
M transmission channels. 

Isolated-resonance approximation 

Owing to the position of the leads' chemical potential relative to the molecular energy levels 
and the large charging energy of small molecules, transport in SMJs is typically dominated by 
individual molecular resonances. In this subsection, we calculate the Green's function in the 
isolated-resonance approximation wherein only a single (non-degenerate or degenerate) molecular 
resonance is considered. In addition to developing intuition and gaining insight into the transport 
mechanisms in a SMJ, we also find (cf. Results and Discussion section) that the isolated-resonance 
approximation can be used to accurately predict the transport. 



N on- degenerate molecular resonance 

If we consider a single non-degenerate molecular resonance then 

r ,^, ,, [n^)+V{^')]C{vy) _ A|A)(A| 

where e = £',y/ — Ey, C{v, v') = A |A)( A| is the rank-1 many-body overlap matrix and we have set 
A = [V^v) + V{i'')]X. In order to solve for G analytically, it is useful to rewrite Dyson's equation 
([6]) as follows: 

G{E) = (1 - G^oiiE)AJ:iE)r' G^,,{E). (10) 

In the elastic cotunneling regime (ASc = 0) we find 

_ A|A)(A[ / A(A|St|A) 

^^^^ - W^TT^ y + E-e + m^ 

A(A|St|A)^^ 



E-e + iQ+ 



A[A)(A| 



^-e-A(A|ST|A) 
Equation [11] can equivalently be expressed as 



(11) 



, in-) + P(--)lg(----). (12) 

where 

t = \V{y)-rV{y')\Tv{C(v,v'yL^] (13) 

is the effective self-energy at the resonance, which includes the effect of many-body correlations 
via the G[y^ v') matrix. 

Using Eq. [H the transmission in the isolated-resonance approximation is given by 

fLfR 

T(E) = ^, (14) 

where 

f ° = [V{u) + V{v')] Tr {C{v, z.')r"} (15) 

is the dressed tunneling-width matrix and F = (F^ + r^)/2. 

As evidenced by Eq. [HI the isolated-resonance approximation gives an intuitive prediction for 
the transport. Specifically, the transmission function is a single Lorentzian resonance centered 



about £ with a half-width at half-maximum of F. The less-intuitive many-body aspect of the 
transport problem is encapsulated in the effective tunneling- width matrices F", where the overlap 
of molecular many-body eigenstates can reduce the elements of these matrices and may strongly 
affect the predicted transport. 

Degenerate molecular resonance 

The generalization of the above results to the case of a degenerate molecular resonance is 
formally straightforward. For an M-fold degenerate molecular resonance 

A ^^ 

The M degenerate eigenvectors of Gmoi ™ay be chosen to diagonalize St on the degenerate subspace 

\{\i\i:^\\j) = 5ijt, (17) 

and Dyson's equation may be solved as before 

G(E).f^M±Z(f^3^(^. (18) 

i^i E — Eyi + £'y; — S/ 

Although St is diagonal in the basis of |A/), F and F need not be separately diagonal. Conse- 
quently, there is no general simple expression for T[E) for the case of a degenerate resonance, but 
T can still be computed using Eq. [TJ 

In this article we focus on transport through Pt-benzene-Pt SMJs where the relevant molecular 
resonances (HOMO or LUMO) are doubly degenerate. Considering the HOMO resonance of 
benzene 

G(E)» CJu^M^^ g'"'-"'). . (19) 

E — ffHOMO — "^uxv' E — enOMO — '^U2v' 

where i^i^2 £ O5 diagonalize St and Oat is the A^-particle ground state. 

PI-ELECTRON EFFECTIVE FIELD THEORY 

In order to model the degrees of freedom most relevant for transport, we have developed an effec- 
tive field theory of interacting vr-electron systems (vr-EFT) as described in detail in Ref.l6. Briefly, 
this was done by starting with the full electronic Hamiltonian of a conjugated organic molecule 
and dropping degrees of freedom far from the vr-electron energy scale. The effective vr-orbitals were 



then assumed to possess azimuthal and inversion symmetry, and the effective Hamiltonian was 
required to satisfy particle-hole symmetry and be explicitly local. Such an effective field theory 
is preferable to semiempirical methods for applications in molecular junctions because it is more 
fundamental, and hence can be readilly generalized to include screening of intramolecular Coulomb 
interactions due to nearby metallic electrodes. 

Effective Hamiltonian 

This allows the effective Hamiltonian for the vr-electrons in gas-phase benzene to be expressed 
as 

"■ {n,m),a 

+ 2YI ^^rnipn - l)(Pm - 1), (20) 

nm 

where t is the tight-binding matrix element, p is the molecular chemical potential, Unm is the 
Coulomb interaction between the electrons on the nth and mth vr-orbitals, and p„ = "^^dnadna- 
The interaction matrix Unm is calculated via a multipole expansion keeping terms up to the 
quadrupole-quadrupole interaction: 

'Jnm ^ UnnUnm 

+ (l-6 ) (u^^^ + u^^ + r/'3*^ + u'^'^) 

where JJ is the monopole-monopole interaction, Jj'^ is the quadrupole-monopole interaction, 
and [/'3'3 is the quadrupole-quadrupole interaction. For two orbitals with arbitrary quadrupole 
moments Qn and Q^ separated by a displacement r, the expressions for these are 

ttMM _ f_ 

er 



where 



ijkl 



Wijki = 5ii5kj + 6ki5ij - 5r ^{rkSurj + rknSij + SkirjVi 
+ri6kjri + Tkridij) + 2>br^*rirjrirk 
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FIG. 1. Spectral functions A{E) = — I/ttTi {G{E)} at room temperature for gas-phase benzene (top panel) 
and Pt-benzene-Pt junctions (ensemble average, bottom panel). The gas-phase resonances are broadened 
artificially as a guide to the eye. The dashed orange lines arc fixed by (left to right) the lowest-l ying op tical 



excitation of the molecular cation 



17 



21| , the vertical ionization energy of the neutral molecule 
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m 



22|, 



and the vertical electron affinity of the neutral molecule [23| . The asymmetry in the average spectral function 
arises because the HOMO resonance couples more strongly on average to the Pt tip atoms than does the 
LUMO resonance. The work function of the Pt(lll) surface (— 5.93eV [241) is shown for reference. 



is a rank-4 tensor that characterizes the interaction of two quadrupoles and e is a dielectric constant 
included to account for the polar izability of the core and a electrons. Altogether, this provides an 
expression for the interaction energy that is correct up to fifth order in the interatomic distance. 

Benzene 



The adjustable parameters in our Hamiltonian for gas-phase benzene are the nearest-neighbor 
tight-binding matrix element t, the on-site repulsion U, the dielectric constant e, and the vr-orbital 
quadrupole moment Q. These were renormalized by fitting to experimental values that should be 
accurately reproduced within a vr-electron only model. In particular, we simultaneously optimized 
the theoretical predictions of 1) the six lowest singlet and triplet excitations of the neutral molecule, 
2) the vertical ionization energy, and 3) the vertical electron affinity. The optimal parametrization 
for the TT-EFT was found to be t = 2.70 eV, U = 9.69 eV, Q = -0.65 eA^ and e = 1.56 with a 



10 



RMS relative error of 4.2 percent in the fit of the excitation spectrum. 

The top panel of Fig. n shows the spectral function for gas-phase benzene within vr-EFT, along 
with experimental values for the first optical excitation of the cation (3.04 eV), the vertical ion- 
ization energy (9.23 eV), and the vertical electron affinity (—1.12 eV). As a guide to the eye, the 
spectrum has been broadened artificially using a tunneling-width matrix of Tnm = (0.2 eY)6nm- 
The close agreement between the experimental values and the maxima of the spectral function 
suggests our model is accurate at this energy scale. In particular, the accuracy of the theoretical 
value for the lowest optical excitation of the cation is noteworthy, as this quantity was not fit 
during the renormalization procedure but rather represents a prediction of our vr-EFT. 

In order to incorporate screening by metallic electrodes into vr-EFT, we utilized an image multi- 
pole method whereby the interaction between an orbital and image orbitals are included up to the 
quadrupole-quadrupole interaction in a screened interaction matrix Unm- In particular, we chose 
a symmetric Unm that ensures the Hamiltonian gives the energy required to assemble the charge 
distribution from infinity with the electrodes maintained at fixed potential, namely 

<Jnm = Unm + Onm^nn + l^i^ ~ "nm){Unm ' ^mn)^ 

(i) 

where Unm is the unscreened interaction matrix and Uiim is the interaction between the nth orbital 
and the image of the rn-th orbital. When multiple electrodes are present, the image of an orbital 
in one electrode produces images in the others, resulting in an effect reminiscent of a hall of 
mirrors. We deal with this by including these "higher order" multipole moments iteratively until 
the difference between successive approximations of Unm drops below a predetermined threshold. 
In the particular case of the Pt-benzene-Pt junction ensemble described in the next section, the 
electrodes of each junction are modeled as perfect spherical conductors. An orbital with monopole 
moment q and quadrupole moment Q^^ located a distance r from the center of an electrode with 
radius R then induces an image distribution at f = — with monopole and quadrupole moments 

and 

5 



^^'■ = -(7) E^^'^^.'Q' 



)kl 
kl 

respectively. Here To. is a transformation matrix representing a reflection about the plane normal 
to the vector r: 

Tik = Sik - 2fifk- 
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FIG. 2. Calculated van der Waals contribution to the binding energy of benzene adsorbed on a Pt(lll) 
surface as a function of distance. Here the plane of the molecule is oriented parallel to the Pt surface. A 
phenomenological short-range repulsion oc r^^'^ has been included to model the Pauli repulsion when the 
TT-orbitals overlap the Pt surface states. 

The lower panel of Fig. „ shows the Pt-benzene-Pt spectral function averaged over the ensemble 
of junctions described in the next section using this method. Comparing the spectrum with the 
gas-phase spectral function shown in the top panel of Fig. „ we see that screening due to the nearby 
Pt tips reduced the HOMO-LUMO gap by 33% on average, from 10.39eV in gas-phase to 6.86eV 
over the junction ensemble. 

The screening of intramolecular Coulomb interactions by nearby conductor(s) illustrated in 
Fig. n leads to an attractive interaction between a molecule and a metal surface (van der Waals 
interaction) . By diagonalizing the molecular Hamiltonian with and without the effects of screening 
included in Unm, it is possible to determine the van der Waals interaction at arbitrary temperature 
between a neutral molecule and a metallic electrode by comparing the expectation values of the 
Hamiltonian in these two cases: 



AK 



vdW 



(H) - {H) 



This procedure was carried out at zero temperature for benzene oriented parallel to the surface of 
a planar Pt electrode at a variety of distances, and the results are shown in Fig. [2l Note that an 
additional phenomenological short-range repulsion oc r~^^ has been included in the calculation to 
model the Pauli repulsion arising when the benzene vr-orbitals overlap the Pt surface states. 

THE LEAD-MOLECULE COUPLING 



When an isolated molecule is connected to electrodes and a molecular junction is formed, the en- 
ergy levels of the molecule are broadened and shifted as a result of the formation of a lead-molecule 
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FIG. 3. The trace of T" for a Pt electrode in contact with a benzene molecule. Nine total basis states of the 
Pt tip are included in this calculation (one s, three p and five d states). The tip height above the plane of the 
molecule is adjusted at each point such that the Pt-C distance is fixed to 2.65A (see text). Tr{r"} retains 
the (six-fold) symmetry of the molecule and is sharply peaked near the center of the benzene ring indicating 
the strongest bonds are formed when the lead is in the 'atop' configuration. The benzene molecule is shown 
schematically with the black lines; the carbons atoms are located at each vertex. 

bond and the electrostatic influence of the leads. The bonding between lead a and the molecule is 
described by the tunneling width matrix T°^ and the electrostatics, including intramolecular screen- 
ing and van der Waals effects, are described by the effective molecular Hamiltonian derived using 
the aforementioned tt-EFT. Although we use the Pt-benzene-Pt junction as an example here, the 
techniques we discuss are applicable to any conjugated organic molecular junction. 

Bonding 



The bonding between the tip of electrode a with the molecule is characterized by the tunneling- 
width matrix F" given by Eq.Sl When a highly-conductive SMJ [J] is formed the lead and molecule 
are in direct contact such that the overlap between the vr-electron system of the molecule and all of 
the atomic-like wavefunctions of the atomically-sharp electrode are relevant. In this case we may 
express the elements of F"^ as 



le{s,p,d,...} 



(21) 



where the sum is over evanescent tunneling modes emanating from the metal tip, labeled by their 
angular momentum quantum numbers, pf{E) is the local density of states on the apex atom of 
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FIG. 4. Eigenvalue decomposition of an ensemble of F" matricies, shovifing that each lead-molecule contact 
has ~ 5 channels. Note that nine orthogonal basis orbitals were included in the calculation for each lead. 



electrode a, and Vp is the tunneling matrix element of orbitals / [25]. The constants C/ can 
in principle be determined by matching the evanescent tip modes to the wavefunctions within 



the metal tip 



251 ]: however, we set Ci = C\ll and determine the constant C by fitting to the 



peak of the experimental conductance histogram [4]. In the calculation of the matrix elements, 
we use the effective Bohr radius of a vr-orbital a* = a^/Z , where oq ~ 0.53A is the Bohr radius 
and Z = 3.22 is the effective hydrogenic charge associated with the vr-orbital quadrupole moment 
-0.65eA2 determined by vr-EFT. 

For each Pt tip, we include one s, three p and five d orbitals in our calculations, which represent 
the evanescent tunneling modes in free space outside the apex atom of the tip. At room tempera- 
ture, the Pt density of states fPOS) p^iE) = ^^ pf{E) is sharply peaked around the Fermi energy 



26l | with p"(ej7')=2.88/eV (27|. In accordance with Ref. l25J, we distribute the total DOS such that 
the s orbital contributes 10%, the p orbitals contribute 10%, and the d orbitals contribute 80%. 

We are interested in investigating transport through stable junctions where the 'atop' binding 
configuration of benzene on Pt has the largest binding energy 28l430[|. In this configuration, the 
distance between the tip atom and the center of the benzene ring is ~ 2.25A [4], giving a tip to 
orbital distance of ~ 2.65A (the C-C bonds are taken as 1.4A). The trace of T°'{eF) is shown as a 
function of tip position in Fig. [3l where for each tip position the height was adjusted such that the 
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distance to the closest carbon atom was 2.65A. From the figure, it is evident that the lead-moIecuIe 
coupHng strength is peaked when the tip is in the vicinity of center of the benzene ring (whose 
outHne is drawn schematically in black). As shown in Ref. [4], the hybridization contribution to 
the binding energy is 

d^Tr{r(^)C(z^,0^)} 






^,/^ y-oo27r -E-E^, + E,, ' 

which is roughly oc Tr{r(ei7)}. Here fi is the chemical potential of the lead metal, T-Ln is the 
A^-particle molecular Hilbert space, and Oat is the ground state of the A^-particle manifold of the 
neutral molecule. The sharply peaked nature of Tr {F"} seen in Fig. [3] is thus consistent with the 
large binding energy of the atop configuration. 

This result motivates our procedure for generating the ensemble of junctions, where we consider 
the tip position in the plane parallel to the benzene ring as a 2-D Gaussian random variable with 
a standard deviation of 0.25A, chosen to corresponded with the preferred bonding observed in this 
region. For each position, the height of each electrode (one placed above the plane and one below) 
is adjusted such that the closest carbon to the apex atom of each electrode is 2.65A. Each lead is 
placed independently of the other. This procedure ensures that the full range of possible, bonded 
junctions are included in the ensemble. 

The eigenvalue distributions of F" over the ensemble are shown in FigUl Although we include 
nine (orthogonal) basis orbitals for each lead, the F matrix only exhibits five nonzero eigenvalues, 
presumably because only five linear combinations can be formed which are directed toward the 
molecule. Although we have shown the distribution for a single lead, the number of transmission 
channels for two leads, where each F" matrix has the same rank, will be the same even though the 
overall lead-molecule coupling strength will be larger. The average coupling per orbital with two 
electrodes is shown in the bottom panel of Fig. [5l 

Screening 

The ensemble of screened interaction matrices Unm is generated using the same procedure 
discussed above. Each Pt electrode is modelled as a conducting sphere with radius equal to the Ft 
polarization radius (1.87A). This is equivalent to the assumption that screening is due mainly to 
the apex atoms of each Pt tip. The screening surface is placed such that it lies one covalent radius 
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FIG. 5. The distribution of charging energy (C/„m) (top panel) and Tr{r} (bottom panel) over the ensemble 
described in the text. Here F = F^ + F^ is the total tunneling-width matrix of the junction. The width 
of the Tr {F} distribution is ^ 4x that of the (C/„m) distribution. The peak of the (Unm) and Tr{F}/6 
distributions are 1.58ey and 1.95ey, respectively, suggesting that transport occurs in an intermediate regime 
where both the particle-like and wave-like character of the charge carriers must be considered. 

away from the nearest carbon atom [Ig ]. 

The average over the interaction matrix elements (Unm) defines the "charging energy" of the 
molecule in the junction [16]. The charging energy {Unm) and per orbital Tr {F} distributions are 
shown in the top and bottom panels of Fig. [H respectively, where two electrodes are used in all 
calculations. As indicated by the figure, the Tr {F} /6 distribution is roughly four times as broad 
as the charging energy distribution. This fact justifies using the ensemble-average Unm matrix for 
transport calculations [2], an approximation which makes the calculation of thousands of junctions 
computationally tractable. The peak values of the (Unm) and Tr{F}/6 distributions are 1.58eV 
and 1.95eV, respectively, suggesting that transport occurs in an intermediate regime where both 
the particle-like and wave-like character of the charge carriers must be considered. 

In addition to sampling various bonding configurations, we also consider the ensemble of junc- 
tions to sample all possible Pt surfaces. The work function of Pt ranges from 5.93eV to 5.12eV for 
the (111) and (331) surfaces, respectively 2j], and we assume that /ipt is distributed uniformly 
over this interval. 

Using this ensemble, the conductance histogram over the ensemble of junctions can be computed, 
and is shown in Fig. [6l The constant prefactor C appearing in the tunneling matrix elements 25 1 
in Eg . [2T] was determined by fitting the peak of the calculated conductance distribution to the that 
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Conductance [Gq] 

FIG. 6. Calculated conductance histogram for the ensemble over bonding configurations and Pt surfaces. 
The value of the conductance peak has been fit to match the experimental data [4] , determining the constant 
C in Eq.[2Il There is no peak for G ^ because we designed an ensemble of junctions where both electrodes 
are strongly bound to the molecule. 

of the experimental conductance histogram [J] . Note that the width of the calculated conductance 
peak is also comparable to that of the experimental peak [J]. 



RESULTS AND DISCUSSION 



The transmission eigenvalue distributions for ensembles of 1.74 x 10^ Pt-benzene-Pt junctions 
calculated using the full many-body spectrum and in the isolated-resonance approximation are 
shown in Figs. [7a] and I7bl respectively. Despite the existence of five covalent bonds between the 

molecule and each lead (cf. Fig.U]), there are only two dominant transmission channels, which arise 

[ i 
from the two- fold degenerate HOMO resonance closest to the Pt Fermi level p]. As proof of this 

point, we calculated the transmission eigenvalue distribution, over the same ensemble, using only 

the HOMO resonance in the isolated-resonance approximation (Eq ll9p . The resulting transmission 

eigenvalue distributions, shown in Fig. I7bl are nearly identical to the full distribution shown in 

Fig. [Taj with the exception of the small but experimentally resolvable J] third transmission channel. 

The lack of a third channel in the isolated-resonance approximation is a direct consequence 

of the two-fold degeneracy of the HOMO resonance, which can therefore contribute at most two 

transmission channels. The third channel thus arises from further off-resonant tunneling. In fact, 

we would argue that the very observation of a third channel in some Pt-benzene-Pt junctions [J] is 

a consequence of the very large lead-molecule coupling (~ 2eV per atomic orbital) in this system. 

Having simulated junctions with electrodes whose DOS at the Fermi level is smaller than that of 

Pt, we expect junctions with Cu or Au electrodes, for example, to exhibit only two measurable 

transmission channels. 
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(a) Full many-body MDE theory (b) Isolated (HOMO) resonance approximation 

FIG. 7. The calculated eigenvalue distributions for an ensemble of 1.74 x 10^ (2000 bonding conhguration 
X 87 Pt surfaces) Pt-benzene-Pt junctions using many-body theory with (a) the full spectrum and (b) the 
isolated-resonance approximation for the (doubly degenerate) HOMO resonance. Despite each lead forming 
~ 5 bonds (cf. Fig. 2]), calculations in both cases exhibit only two dominant channels which arise from the 
degeneracy of the relevant (HOMO) resonance. The weak third channel seen in (a) is a consequence of the 
large lead-molecule coupling and is consistent with the measurements of Reflj. 

In order to investigate the efficacy of the isolated-resonance approximation further, we calcu- 
lated the average total transmission through a Pt-benzene-Pt junction. The transmission spectra 
calculated using the full molecular spectrum, the isolated HOMO resonance and the isolated 
LUMO resonance are each shown as a function of the leads' chemical potential /Upt in Fig. [HI 
The spectra are averaged over 2000 bonding configurations and the blue shaded area indicates 
the range of possible chemical potentials for the Pt electrodes. The close correspondence between 
the full transmission spectrum and the isolated HOMO resonance over this range is consistent 
with the accuracy of the approximate method shown in Fig. [71 Similarly, in the vicinity of the 
LUMO resonance, the isolated LUMO resonance approximation accurately characterizes the av- 
erage transmission. The HOMO-LUMO asymmetry in the average transmission function arises 
because the HOMO resonance couples more strongly on average to the Pt tip atoms than does 
the LUMO resonance. 



It is tempting to assume, based on the accuracy of the isolated-resonance approximation in our 
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FIG. 8. The calculated average total transmission averaged over 2000 bonding configurations through a 
Pt-benzene-Pt junction shown as a function of the leads' chemical potential /ipt- The isolated resonance 
approximation using the HOMO or LUMO resonance accurately describe the full many-body transport in 
the vicinity of the HOMO or LUMO resonance, respectively. These data are in good agreement with the 
measurements of Ref. Ij. The work-function range for the crystal planes of Pt are shaded in blue, where 
-5.84ey < fj.pt < -5.12ey [24|. 



many-body transport theory, that an analogous "single molecular orbital" approximation would 
also be sufficient in a transport calculation based e.g. on density-functional theory (DFT). However, 
this is not the case. Although the isolated-resonance approximation can also be derived within 
DFT, in practice, it is necessary to use an "extended molecule" to account for charge transfer 
between molecule and electrodes. This is because current implementations of DFT fail to account 
for the particle aspect of the electron 



10|. Analyzing transport in terms of extended molecu 
orbitals has proven problematic. For example, the resonances of the extended molecule in Ref. 
apparently accounted for less than 9% of the current through the junction. 



Using an "extended molecule" also makes it difficult, if not impossible, to interpret transport 



contributions in terms of the resonances of the molecule itself [31[. Since charging effects in SMJs are 
well-described in our many-body theory g, II4], there is no need to utilize an "extended molecule," 
so the resonances in our isolated-resonance approximation are true molecular resonances. 



The full counting statistics of a distribution are characterized by its cumulants. Using a single- 



particle theory to describe a single-channel junction, it can be shown 



M 



33t | that the first cumulant 



is related to the junction transmission function while the second cumulant is related to the shot 
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FIG. 9. The calculated Fano factor F distribution for the full ensemble of 1.74 x 10^ Pt-benzene-Pt 
junctions. F describes the nature of the transport, where F = {) and F =\ characterize wave-like (ballistic) 
and particle-like transport, respectively. The peak value of this distribution F ^ 0.51 indicates that we are 
in an intermediate regime. 



(22) 



noise suppression. Often this suppression is phrased in terms of the Fano factor [3j] 

p Y^ Tn{^ — Tn) 

In Fig. [9] we show the distribution of F for our ensemble of junctions, where the r„ have been 
calculated using many-body theory. Because of the fermionic character of the charge carriers 
< -F < 1, with F = corresponding to completely wave- like transport and a value oi F = \ 
corresponding to completely particle-like transport. From the figure, we see that F is peaked 
~ 0.51 implying that the both particle and wave aspects of the carriers are important, a fact which 
is consistent with the commensurate charging energy and bonding strength (cf. Fig. [5|). 

In such an intermediate regime both the 'complementary' aspects of the charge carriers are 
equally important, requiring a many-body description and resulting in many subtle and interesting 
effects. For example, the transport in this regime displays a variety of features stemming from the 
interplay between Coulomb blockade and coherent interference effects, which occur simultaneously 
[O, 111]. Although the Fano factor reflects the nature of the transport, it is not directly related to 
the shot-noise power in a many-body theory. The richness of the transport in this regime, however, 
suggests that a full many-body calculation of a higher-order moment, such as the shot-noise, may 
exhibit equally interesting phenomena. 



CONCLUSION 



We have developed a state-of-the-art technique to model the lead-molecule coupling in highly- 
conductive molecular junctions. The bonding between the lead and molecule was described using an 
'ab initio' model in which the tunneling matrix elements between all relevant lead tip wavefunctions 
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and the molecule were included, producing multi-channel junctions naturally from a physically 
motivated ensemble over contact geometries. Coulomb interactions between the molecule and the 
metallic leads were included using an image multipole method within vr-EFT. In concert, these 
techniques allowed us to accurately model SMJs within our many-body theory. 

The transport for an ensemble of Pt-benzene-Pt junctions, calculated using our many-body 
theory, confirmed our statement [2] that the number of dominant transmission channels in a SMJ 
is equal to the degeneracy of the molecular orbital closest to the metal Fermi level. We find 
that the transport through a Pt-benzene~Pt junction can be accurately described using only the 
relevant (HOMO) molecular resonance. The exceptional accuracy of such an isolated-resonance 
approximation, however, may be limited to small molecules with large charging energies. In larger 
molecules, where the charging energy is smaller, further off-resonant transmission channels are 
expected to become more important. 

In metallic point contacts the number of channels is completely determined by the valence of the 
metal. Despite the larger number of states available for tunneling transport in SMJs, we predict 
that the number of transmission channels is typically more limited than in single-atom contacts 
because molecules are less symmetrical than atoms. Channel-resolved transport measurements of 
SMJs therefore offer a unique probe into the symmetry of the molecular species involved. 
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